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^r^ Abstract — We have developed a method for analyzing the kinematic association of isolated relativistic 

^""l objects— possible remnants of disrupted close binary systems. We investigate pairs of fairly young radio 

,__, pulsars with known proper motions and estimated distances (dispersion measures) that are spaced no more 

QCi, than 2-3 kpc apart. Using a specified radial velocity distribution for these objects, we have constructed 

C/^ 100-300 thousand trajectories of their possible motion in the Galactic gravitational field on a time scale of 

,^ several million years. The probabilities of their close encounters at epochs consistent with the age of the 

»p^ younger pulsar in the pair are analyzed. When these probabilities exceed considerably their reference values 

O obtained by assuming a purely randomencounter between the pulsars under consideration, we conclude that 

-4— > the objects may have been gravitationally bound in the past. As a result, we have detected six pulsar pairs 

^ (J0543+2329/J0528+2200, J1453-6413/J1430-6623, J2354+6155/J2321+6024, J1915+1009/J1909+1102, 

' ' J1832-0827/J1836-1008, and J1917+1353/J1926+1648) that are companions in disrupted binary systems 

T-H with a high probability. Estimates of their kinematic ages and velocities at binary disruption and at the 

^ present epoch are provided. 
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INTRODUCTION 

Neutron stars (NSs) are among the fastest Galactic 
objects. Their tangential velocities determined from 
measured proper motions reach hundreds or even 
thousands of kilometers per second (see, e.g., Hobbs 
et al. (2005) and references therein). NSs probably 
receive such kicks during asymmetric supernova 
explosions (Iben and Tutukov 1996; Helfand and 
Tademaru 1977; Shklovsky 1970; Dewey and Cordes 
1987) or due to disruption by the explosions of 
the binary systems whose components they were 
(Blaauw 1961; Gott et al. 1970). The third possibility 
is the birth of NSs during the explosions of 
highvelocity massive stars that have been ejected 
from young star clusters by dynamical processes in 
their cores (Poveda et al. 1967; Gvaramadze 2007; 
Gvaramadze et al. 2008, 2009). 

At least 50% of all stars are known to be 
members of binary and multiple systems (Batten 
1967; Zasov and Postnov 2006; Duquennoy and 
Mayor 1991; Halbwachs et al. 2003) and, accordingly, 
many of the relatively young NSs that we observe as 
pulsars must have been formed precisely in binary 
systems. 

NSs are the final evolutionary stage of stars 
with initial masses in the range from 8-10 to 25-30 
Mq (see, e.g., Postnov and Yungelson 2006; Bethe 
and Brown 1998), although there is evidence that 
the initial masses can be even higher (Muno et al. 
2006). Most such close binary systems (CBSs) are 
disrupted by a supernova explosion, for which the 
loss of half of itsmass is sufficient in the simplest case 
of a symmetric explosion and a circular orbit; for 
an asymmetric explosion, this fraction can be even 
smaller. 

Nevertheless, when one of the stars in a 
massive binary system explodes as a supernova at 
the end of its life, this does not necessarily lead 
to binary disruption — about 40% of the pairs 
withstand the explosion (Bethe and Brown 1998). 
Since the components of a massive CBS exchange 
mass during its evolution, only the C-0 core might 
remain of the primary (more massive) star at the 
time of its collapse. For this reason, the mass lost 



during the first supernova explosion may prove to be 
insufficient for binary disruption; only the orbits of 
the components become eccentric and the center of 
mass acquires a velocity of about 100 km s~^. A CBS 
with a relativistic component, a NS or a black hole 
(BH), emerges (see, e.g., Zasov and Postnov 2006). 

When the secondary component collapses 
during the subsequent evolution (the second 
supernova explosion occurs) and the second 
relativistic object is formed, the pair is generally 
disrupted: the binary as a whole is already less 
massive and, hence, the ejected mass (higher than 
that in the first case) most often exceeds the value 
needed for disruption (Zasov and Postnov 2006). 

The high velocities of pulsars serve as an 
additional argument that many of them were 
members of disrupted binary systems in the past. 
The slingshot effect acts during the disruption of a 
binary system, and this is one of the most natural 
mechanisms for an increase in the space velocities 
of disruption remnants. The kick received by the 
pulsar during an asymmetric supernova explosion in 
conjunction with the orbital angular momentum can 
accelerate the NS significantly (Bailes 1989; Kiel and 
Hurley 2009). However, this mechanism does not act 
at every CBS disruption. 

These currently isolated NSs may prove 
to be a convenient tool for investigating the 
disruptionmechanism. However, the kinematic 
characteristics of isolated pulsars per se do not allow 
even the contribution from additional velocities of a 
different nature to be estimated. The situation will 
change radically if we were to find two pulsars that 
used to be companions in a disrupted CBS. 

Considering pairs of fairly close pulsars with 
known proper motions and distances and specifying 
the distribution of their space velocities, we can trace 
themotion of the objects in the past and determine 
the distance of their closest approach. If its epoch 
turns out to be consistent with the characteristic age 
Teh of one of the pulsars {tch = P/2P, where P is the 
pulsar period) , then this pair can be assumed to have 
a common origin. The significance of this hypothesis 
can be estimated by comparing the probability of 



the encounter found with that for any two objects 
with arbitrary proper motions. The presence of a 
supernova remnant in the corresponding Galactic 
region that is also consistent in age with one of 
the pulsars can serve as an additional criterion 
for the validity of the identification. Accepting the 
hypothesis about a common origin will allow us to 
analyze the collapse and second supernova explosion 
kinematics and the CBS disruption dynamics as 
well as to ascertain the relationship between the 
characteristic and actual pulsar ages. 

Obviously, it is possible to establish 
membership in the same CBS only at disruption 
in the second supernova explosion. The epoch of 
closest approach must then coincide with the age 
of the younger pulsar. If, however, the binary was 
disrupted after the first explosion, then the second 
(intact) star evolves independently and can recede to 
hundreds of parsecs before it collapses 10^ — 10^ yr 
later. The pulsar born in this asymmetric explosion 
that received an additional kick will change the 
direction of its motion, which will make the detection 
of its companion virtually impossible. Thus, the 
total number of disrupted pairs can be larger. 

Another scenario inaccessible to analysis by 
the proposed method is that the massive progenitor 
stars of pulsars in the past were ejected from young 
clusters and, accordingly, acquired a high space 
velocity (Gvaramadze 2007; Gvaramadze et al. 2008, 
2009). In that case, they evolve and collapse as 
isolated stars. If the explosion does not change 
greatly the direction of their velocities, then tracing 
the trajectories of the two emerged young NSs into 
the past will only point to the same cluster, i.e., the 
encounter will be fairly close, but the stars cannot 
be said to have been gravitationally bound in the 
past. Therefore, on the one hand, the pairs of pulsars 
whose trajectories start near young star clusters 
should be dealt with carefully (see, e.g., Vlemmings 
et al. 2004) and, on the other hand, it should be 
understood that the described situation is unlikely. 
The progenitor stars ejected from clusters must be 
massive enough to subsequently explode as SNe II. 
Only about 20 fast massive stars are known to date, 
and only two of them have a mass higher than 8 Mq 



(Przybilla et al. 2008; Heber et al. 2008). 

The age difference between the components 
(i.e., the time elapsed between the first and 
second supernova explosions) for massive binaries 
should not exceed several million years. The latter 
requirement follows from the characteristic lifetime 
of massive main-sequence stars. If the time between 
the explosions is short, then the characteristic ages 
of both components can coincide. 

The question about the actual ages of pulsars 
remains an open one. The characteristic age, 
particularly for young pulsars, can serve only as 
a reference. Faucher-Giguere and Kaspi (2006) 
analyzed the relationship between the actual and 
characteristic ages and showed that they could differ 
by a factor of several. Additionally, Teh does not need 
to be an upper limit for the actual age (see, e.g., 
Gaensler and Frail 2000). The described approach to 
identifying companion pulsars will help to estimate 
their ages more accurately. 

Isolated BHs that were previously members of 
disrupted CBSs together with pulsars can also be 
searched for by a similar method. The BH is the 
first to be formed at the end of the evolution of 
the more massive (> 30 Mq) star, following which 
the binary survives with a probability of 40% and is 
disrupted after the second supernova explosionwith 
the formation of a NS (Bethe and Brown 1998). In 
this case, since the BH is massive, it will acquire 
a lower velocity and, hence, will be not far from 
the place of CBS disruption (the pulsar birthplace). 
Given the kinematic characteristics of the pulsar (as 
one of the components of the disrupted binary), 
the BH localization region where it makes sense 
to search for it as a peculiar source of radiation 
with various frequencies can be limited in space 
(and in the sky) (Beskin and Karpov 2005). The 
remnant of the second supernova near which the 
pulsar trajectory passes if their ages are consistent, 
can be an important reference in this search. 

To localize the BH in such a way, we suggest 
using pulsars with accurately measured proper 
motions and distances. This accuracy will directly 
determine the size of the region (from 1 to 10 



square degrees) where the BH is located under 
the assumption that it exists. In these regions, 
the catalogs of starlike sources in all ranges will 
be cross identified and peculiar objects (primarily 
those without any lines in their optical spectra) 
will be selected to search for superfast variability 
in them, a basic property of an isolated accreting 
BH (Beskin et al. 1997; Beskin and Karpov 2005). 
Pulsars with young characteristic ages will be used, 
where possible. In this case, the spatial separation 
of the objects (the BH and the pulsar) is definitely 
small (even given the inaccuracy of the age estimate) , 
which will facilitate the search for the BH. 

Note separately that runaway OB stars can 
also be used to search for NSs and BHs (see 
Prokhorov and Popov (2002) and references therein). 

In this work, we search for the pairs of isolated 
radio pulsars that were bound in the past. The 
detection of such pairs will allow one to understand 
the CBS disruption dynamics and to estimate the 
velocities acquired by the pulsars during a supernova 
explosion, their radial velocities, ages, and type of 
explosion. 

Several dozen pulsars associated with 
supernova remnants are known to date. The 
most comprehensive overview of such associations 
can be found in Tian and Leahy (2004). These are 
mostly young pulsars located either within or in the 
immediate vicinity of their remnants. However, we 
are interested in the cases where the pulsar is fast 
and old enough to have left the neighborhood of the 
remnant (observed until now). 

The attempts to discover disrupted pulsar 
pairs were previously made by Bisnovatyi-Kogan and 
Romberg (1974), Gott et al. (1970), and Wright 
(1979). The next step, an analysis of the possible 
association of a pair of moderately young pulsars 
spaced well apart for which the parallaxes and 
proper motions are known with a good accuracy, was 
undertaken by Vlemmings et al. (2004). Below, we 
discuss this result. 

In the next section, we describe in detail the 
methodology and the assumptions made here. Next, 
we present the results of our investigation of specific 



pulsar pairs. The work is completed by a discussion 
of the results obtained. 



THE METHOD 

Formulation of the Problem 

Searching for the kinematic association of two 
pulsars with known proper motions consists 
in determining the closest approaches of their 
reconstructed trajectories in the Galaxy in the past. 
Obviously, the pulsars with a common origin (i.e., 
those that were the components of the same binary 
system) were closest to each other at the time of CBS 
disruption. This epoch corresponds to the current 
age of the younger pulsar (recall that we analyze the 
case where the system was disrupted after the second 
supernova explosion). 

We consider all motions in a rectangular 
coordinate system whose origin is placed at the 
Galactic center. One of its axes is parallel to the 
direction of the Sun, the second axis is directed along 
the velocity of the local standard of rest, and the 
third axis is perpendicular to the Galactic plane and 
complements the first two to a right-handed triple 
of vectors. In such a coordinate system, the radius 
vector of the Sun is Tq = {—8.5 kpc ,0,0} (see, e.g., 
Mihalas and Binney 1981). 

The law of pulsar motion r(t) in the Galactic 
gravitational potential v'G(r) is the solution of the 
equation of motion 

r = -Vipair) (1) 

with the initial conditions 



ro = r(t = 0),Vo = V(t = 0), 



(2) 



where the present epoch corresponds to the initial 
time. For known functions ri^2(i) for the pulsar pair, 
the separation p between them at the past time 
—t is |ri(— t) — r2(— 1)|. Associating the pulsars, we 
seek for such an epoch T consistent with the age 



of the younger pulsar (no older than its several \^r,rot attributable to the Galactic rotation of the 
characteristic ages) that pulsar LSR (i.e., also to the Galactic rotation curve) 

and a peculiar velocity V^^p. Thus, we finally obtain 



|ri(-T) - r2(-r)| = min/9. 



(3) 



Vo = ^fr,rot + V^,p + Vi + rg 



0- 



(7) 



A low value of minp can indicate that these 
objects were gravitationally bound in the past. The Models Used 



Solution of the Equation of Motion 

Since the Galactic gravitational potential is not 
spherically symmetric (see below), the solution of 
Eq. (1) is generally obtained numerically, while the 
vectors tq and Vq are determined directly from 
observations. For the current radius vector of a pulsar 
with Galactic coordinates / and h and a heliocentric 
distance d, we have 



ro = (i{cos6cos/,cos6sin/, sin 6} + yq. (4) 

Its current velocity vector Vq is specified by the 
known proper motion {^i cos h, //{,}, distance d, radial 
velocity T^, and solar velocity vector 



r© 



Vo .ot + Vp 



Q,LSR- 



(5) 



Here, \^Q^rot = {0, Vgi^j-ot) 0} is the Galactic rotation 
velocity of the local standard of rest (LSR). It is 
determined by the Galactic rotation curve (i.e., by 
the potential ipcir)) and is about 220 km s^^. The 
vector ^Q^LSR is the solar velocity relative to the 
LSR; it is directed to the point with I = 53° and 
b = 25° (apex) and its magnitude is 16.5 km s^^ 
(Mihalas and Binney 1981). 

The proper motion of a pulsar specifies its 
transverse velocity V^ whose magnitude (in km s"-*^) 
is 



Vt = 4.74 • d • V(wcos6)2 + 



/Ufe^ 



(6) 



where fii and /if, are in [mas yr ^] and d is in kpc. The 
pulsar radial velocity consists of a secular component 



Apart from the spherical coordinates, only the proper 
motions of pulsars and, more rarely, their distances, 
when the parallax vr has been measured (the 
distances are usually calculated from the dispersion 
measure), are generally accessible to us directly from 
observations. All of these quantities are known with 
a finite accuracy and the distances are most often 
mo del- dependent . 

The main problem in studying the kinematics 
of pulsars is the absence of direct measurements of 
their radial velocities. At the same time, the source 
of information about the latter is the distribution 
function of one (any) of the three peculiar velocity 
components of a pulsar at the time of its birth 
(e.g., at the time of binary disruption). The total 
peculiar velocity is assumed to be isotropic relative 
to the LSR. Therefore, the same function can also be 
used for the peculiar component of the pulsar radial 
velocity. 

It is constructed from the set of data on 
all of the investigated pulsars as a solution of 
the maximum likelihood problem. At present, the 
question about its unambiguous choice is sill an 
open one. Various versions have been proposed: 
Gaussian and non-Gaussian, one- or two-component 
distributions (see, e.g., Hobbs et al. 2005; Cordes and 
Chernoff 1998). Faucher-Giguere and Kaspi (2006) 
compare in detail the results available to date and 
conclude that it is so far hard to unambiguously 
choose a particular model, because all models are 
consistent with the observational data. 

Thus, in the modeling task (in contrast to 
theoretical constructions), the choice of a particular 
distribution is not so important .We used it in the 
form of a sum of two Gaussians (Arzoumanian et al. 
2002): 



p(yr,p) 



Wl 



2iTai 



exp 



+ 



1 — Wl 



exp 




where ai 
0.4. 



^27rcr2 
90 km s~^, (T2 = 500 km s^^ and wi 



(8) 



As was said above, it describes the velocities 
of the pulsars at the epoch of their birth and can 
change by now. It is easy to estimate the variations 
AV in velocity magnitudes as a result of the motion 
in the Galactic potential: 



AV 



r^ rr, 



(9) 



where r is the pulsar age and f is the magnitude 
of its acceleration. Vlemmings et al. (2004) used the 
characteristic age Teh as r and the acceleration z in a 
direction perpendicular to the Galactic plane instead 
of r. However, they investigated a pair of moderately 
young pulsars with ages TcH ~ 3x 10^ yr. In this time, 
correction (9) may prove to be significant and, at the 
same time, significantly different from the true one, 
since Teh is not a good age estimate and because the 
corresponding integral is replaced with product (9). 

Although here we investigate pulsars with 
considerably younger ages, to estimate the effect of 
the Galactic gravitational potential on their motion, 
we simulated sets of pulsar trajectories that began 
their motion in the Galactic plane (where the effect 
will be at a maximum), at Galactocentric distances 
of 2.5 — 10.5 kpc (with a step of 2 kpc). We 
obtained changes in peculiar velocity magnitudes and 
deviations of the trajectories from rectilinearity after 
10^ 5xl0^ 10^ and 3x10^ yr. 

Function (8) was used as the distribution 
of pulsar initial peculiar velocity components. The 
gravitational potential fci^) was specified by an 
expression from Carlberg and Innanen (1987) and 
Kuijken and Gilmore (1989): 



<^Gir,z) 



GMdh 



flG + Y.UP^4^^Y + hlh + r 



Table 1. Parameters of the Galactic potential (Carlberg 
and Innanen 1987; Kuijken and Gilmore 1989) 
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(10) 



where r is the distance from the Galactic center 
and z is the distance from the Galactic plane. This 
is a three-component axisymmetric potential that 
includes the contribution from the disk-halo, bulge, 
and nucleus; its parameters are presented in Table 1. 

Obviously, the closer the pulsar to the Galactic 
center, the more nonuniform and nonrectilinear 
its motion. The limiting case for a Galactocentric 
distance of 2.5 kpc is presented in Figs. 1 and 
2, which shows the changes in peculiar velocity 
magnitudes for pulsars and the deviations of their 
trajectories from a rectilinear shape. We see that 
their motion is almost uniform (compared to 
the parameters of distribution (8)) and rectilinear 
on time scales that do not exceed 10^ yr. All 
distributions are bimodal, because function (8) is a 
two-component one. 

For pulsars that began their motion from 
a Galactocentric distance of 10.5 kpc, the scatter 
0"[Ay] does not exceed 4 km s^-*^ even for T = 3x 10^ 
yr and the mean ;u[Ar] = 21 pc. 

Thus, when the motion of pulsars in the past 
is investigated on time scales of ~ 10^ yr or shorter, 
including the effect of v^clr) in the form of correction 
(9) seems superfluous and can even lead to additional 
errors in the results, since the actual age r can 



differ from the characteristic one Teh by a factor of 
several (see, e.g., Faucher-Giguere and Kaspi 2006). 
In contrast, distribution (8) also describes well the 
current pulsar velocities. 

The initial conditions Tq and Vq are specified 
in the form of distributions that reflect the 
inaccuracies in the distances and velocities. 



Randomization of Initial Conditions 

As has been noted above, all of the pulsar 
characteristics known from observations are specified 
with a finite accuracy. As in Vlemmings et al. 
(2004), we assume the pulsar proper motions to 



w^ 



be distributed according to normal laws N(fii,a 
and N{fih,a'f^) with means equal to the measured 
values and variances equal to the squares of 
themeasurement errors. A similar assumption is 
also used for the distances (mean d and variance 
(T^). Here, we estimated them from highly accurate 
measurements of the dispersion measure (DM) based 
on the NE2001 model of the Galactic free electron 
density distribution (Cordes and Lazio 2002) with a 
relative error of about 30%. 

The pulsar equatorial coordinates are known 
with very high accuracies of fractions of an 
arcsecond, which gives an uncertainty in the position 
of 10~^ — 10~^ pc at distances to the objects of several 
kpc. These values are considerably lower than the 
distance errors aa and, hence, we assumed a and 
S (as well as the calculated / and b) to be known 
exactly. 

Since, as was shown above, the pulsar motion 
on time scales of 10^ yr remains almost uniform, the 
uncertainty in the calculated r(t) will generally be 
the same as that for the distance d. 

In turn, 0"^^ and 0"^,^, along with the error in 
the distance d, will lead to randomization of the 
magnitude (and direction) of the current velocity Vq. 
For d ~ 2 — 4 kpc, proper motions of ~ 10 mas yr~^, 
and their relative errors of ~ 30%, the error in the 
transverse velocity will ~ 100 km s~^, which will 
subsequently, on time scales up to 10^ yr, lead to 
an error in the position r{t) of 100 pc. This is much 



smaller than the contribution from the inaccuracy in 
the initial distance d. 



Simulations 

It is clear from the aforesaid how the distributions 
Pr(ro) and Pvi^o) of initial vectors fq and Vq 
in (2) are obtained. Searching for the closest 
encounters of two pulsars is reduced to constructing 
the distribution of the probability P{p) that the 
minimum separation between them (on a given time 
scale in the past) does not exceed p. 

This distribution was constructed by the 
Monte Carlo method; solving Eq. (1) and simulating 
the sets of pairs of pulsar trajectories with various 
realizations of fq and Vq from Pr(ro) and p«(Vo), 
we traced them into the past to a time of 1-2 
Myr (exceeding the characteristic age of the younger 
pulsar) and calculated the minimum separation pmin 
between them. Since the effect of the Galactic 
gravitational potential (10) in our problem is 
negligible, we restricted ourselves to several points 
with a step of 500 000 yr determined by the 
Runge-Kutta-Feldberg method for each trajectory. 
The pulsar motion between them was assumed to be 
uniform and rectilinear, which allowed the speed of 
our simulations to be increased significantly. 

Together with the distribution P{p), we 
obtained the distribution T{p) for the epochs of 
closest approaches. Next, after comparison of the 
former distribution with the reference one Po{p), we 
reached the conclusion about whether the kinematic 
association of the pulsar pair under study was 
possible. 

Reference Probability 

Once the distribution P{p) has been obtained, we 
need a criterion that would allow us to determine 
how significant the probability that a given pair of 
pulsars approached to some small distance p in the 
past is. If the probability P{p) is too high for the 
simulated number of approaches to a distance ^ p to 
be explained only by random factors, then the 
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pair under study proves to be kinematically singled 
out among the others. There is reason to beheve 
that these pulsars were once actually close to each 
other and to suggest that they constituted a single 
disrupted binary system. 

This necessary criterion consists in comparing 
the probabilities P{p) and Po{p), where the latter is 
the probability for two randomly chosen pulsars to be 
at a distance no larger than p during their motion. 
Obviously, this probability should be calculated by 
taking into account the pulsar density in the Galaxy. 

Since there are pulsars of different ages and 
different origins in the Galaxy, Po{p) is equivalent to 
the probability of detecting two (or more) pulsars 
at a distance ^ p when a given region of space 
is observed. Clearly, the larger this region and the 
higher the pulsar density in it, the higher the 
probability of such a detection. 

Thus, to calculate Po{p), we must specify, first, 
the pulsar density distribution over the Galaxy and, 
second, the sizes and localization of the region of 
space of interest. Here, we used the radial pulsar 
surface density distribution derived by Yusifov and 
Kiiciik (2004): 



S(i?) = A 



exp 



R + Ri 
Rq + i?i 

R — Rq 
Rq + Ri 



:ii) 



where A = 37.6 ± 1.9 kpc^^ a = 1.64 ± 0.11, 
b = 4.01 ± 0.24, Ri = 0.55 ± 0.10 kpc, and Rq 
= 8.5 kpc, and a similar distribution in a direction 
perpendicular to the Galactic plane, whose choice 
and the optimality of the parameter {z) = 50 pc 
are justified by Faucher-Giguere and Kaspi (2006): 



p{z) 



1 



2{z) 



exp 



\z\ 



(12) 



Collectively, these distributions specify the 
average number of pulsars r][r) that can actually be 
detected in a Galactic region with coordinates of its 
center r = {x, y, z} and sizes p [pc]: 



r/(r) = 10-V5^(Va;2+y2) / p{z)dz. (13) 

J z-ol2 



Here, x, y, z are the rectangular coordinates in 
the coordinate system associated with the Galaxy 
and described above. The distribution of pulsars is 
assumed to be stable for several million years. 

Calculating the probability P{p) from our 
simulations, we automatically obtain the region (rp, 
Tp + drp) where the pulsars approach each other 
to distances no larger than p. When p <C drp, the 
number of pulsars in the volume p^ obeys a Poisson 
distribution. Therefore, we divided the entire region 
into nonoverlapping cells r^ (A; = l..A^) with sizes p 
in each of which the probability of detecting fewer 
than two pulsars is specified by the expression 



e^ = E 



Wk)]' 



n<2 



n\ 



exp[-r?(rfc)]. 



(14) 



Finally, Pq (p) ) the probability of detecting two 
or more pulsars at a distance ^ p somewhere in the 



Galactic region 
the formula 



^p' 



Tp + dYp) 



was calculated from 



N 



p«m 



fc=i 



(15) 



This quantity is indicative of whether a 
particular pulsar pair stands out kinematically. For 
many pairs, P{p) and Po{p) were comparable (see 
below), but they differed significantly for some of 
them. In our problem, the kinematic association of 
pulsars is based on Monte Carlo calculations of the 
probability that a specific pair of pulsars approached 
to a given (or smaller) distance in the past. If this 
probability for sufficiently close approaches (of the 
order of several parsecs) is high compared to the 
reference one P{p) ^ Poip)i then it is concluded 
that they have a common origin. 




500 



Fig. 3: Distribution of minimum separations between the 
pulsars J1453-6413 and J1430-6623 derived from our 
simulations of 15 000 trajectories. 



Mutual Correspondence 
Adopted Models 



between the 



Here, we assume the distance to the Galactic center 
to be Rq — 8.5 kpc. The gravitational potential 
(10) suggests the same Rq. It is based on a 
standard Galactic rotation curve. In contrast, the 
distribution of pulsar peculiar velocities obtained by 
Arzoumanian et al. (2002) uses the TC93 models 
of the Galactic free electron density (Taylor and 
Cordes 1993), while we use NE2001. Their spatial 
distribution is assumed to be uniform and the 
gravitational potential is based on the papers by 
Paczynski (1990) and Caldwell and Ostriker (1981), 
where Rq — 8 kpc. In general, as is pointed out in 
the "Models Used" Section, the observational data 
available to date do not allow one to unambiguously 
choose a particular distribution of pulsar velocities. 
We are talking not only about the specific parameters 
but also about the shape of the distribution (for 
a review, see Faucher-Giguere and Kaspi (2006)). 
All of the possible versions do not differ from one 
another dramatically. For our results to change 
qualitatively, it is necessary that the simulated 
probabilities of close encounters be several orders of 
magnitude higher than the values obtained (as will 
be seen from the results). Varying the shape of the 
velocity distribution cannot lead to such a situation. 



especially since the pulsar motion is investigated on 
time scales of several hundred thousand years, when 
reasonable differences in the distributions do not yet 
manifest themselves. Therefore, the distribution of 
peculiar velocities from Arzoumanian et al. (2002) is 
quite suitable for this work. 



RESULTS 

Calculations 

The possibility of a kinematic association of pulsars 
was considered for pairs that satisfied the following 
conditions: 

(1) the proper motions (^a, w) ^^^d distances 
(parallax vr, DM) are known; 

(2) the current separation does not exceed 2-3 
kpc; 

(3) the characteristic age T^h = P/2P of the 
younger pulsar does not exceed 10^ yr (this is because 
the binary disruption region is determined more 
accurately for a small displacement of the pulsar from 
its current position backward along the trajectory, 
especially since we can be in error by a factor of 
several while estimating the disruption time from the 
characteristic age); 

(4) the pulsar age difference does not 
exceed several million years (this follows from 
the characteristic lifetime of massive main-sequence 
stars). If the true pulsar age is no more than ~100 
000 yr, then the supernova remnant associated with 
it can survive (and be detected) to the present day. 

Here, we initially formed a sample of ~200 
isolated radio pulsars with known proper motions 
and distances. Their parameters were taken mostly 
from the ATNF database (Manchester et al. 2005; 
http://www.atnf.csiro .au/research/pulsar/psrcat/expert.html) 
We selected 16 pairs of objects from them that met 
the above criteria. For each of these pairs, we first 
simulated 10-15 thousand trajectories (per each 
pulsar). Next, we selected 7 from the 16 pairs for 
which the probability density p{p) = -^P{p) of the 
distribution of separations between the components 
had a distinct sharp peak at p ^ 500 pc. 
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Table 2. Pulsar parameters. 



Pulsars 


a 
(lims) 


5 
(dms) 


D 

(kpc) 


(mas yr"^) 


(mas yr^-*^) 


Teh 

(10^ yr) 


A 
(kpc) 


J0543+2329 
J0528+2200 


05 43 09.660 
05 28 52.308 


23 29 05.00 
22 00 01.00 


2.06 
1.61 


19(7) 
-20(19) 


12(8) 
7(9) 


253 

1480 


0.46 


J1453-6413 
J1430-6623 


14 53 32.737 
14 30 40.872 


-64 13 15.59 
-66 23 05.04 


2.08 
1.00 


-16(1) 
-31(5) 


-21.3(8) 
-21(3) 


1040 
4490 


1.08 


J2354+6155 

J2321+6024 


23 54 04.724 
23 21 55.213 


61 55 46.79 
60 24 30.71 


3.43 
3.03 


22(3) 
-17(22) 


6(2) 
-7(19) 


920 
5080 


0.46 


J1915+1009 
J1909+1102 


19 15 29.98290 
19 09 48.69380 


10 09 43.780 

11 02 03.350 


6.27 
4.18 


4(7) 
-6(4) 


-10(13) 

7(8) 


420 
1700 


2.09 


J1832-0827 
J1836-1008 


18 32 37.0200 
18 36 53.925 


-08 27 03.64 
-10 08 08.3 


4.85 
4.46 


-4(4) 
18(65) 


20(15) 
12(220) 


161 
756 


0.42 


J1833-0827 
J1836-1008 


18 33 40.3000 
18 36 53.925 


-08 27 31.25 
-10 08 08.3 


4.66 
4.46 


-38(26) 
18(65) 


-72(55) 
12(220) 


147 
756 


0.24 


J1917+1353 
J1926+1648 


19 17 39.7902 
19 26 45.322 


13 53 56.95 
16 48 32.78 


3.99 

5.83 


0(12) 
13(11) 


-6(15) 
-14(18) 


428 
511 


1.87 



Table 3. Pairs of pulsars that approached each other most closely in the past. 



Pulsars 


N 


Po 
(pc) 


Pmin 

(pc) 


n(5) 


P(5) 


^o(5) 


PoiPo) 


T(5) 
(103 yr) 


K-(5) 
(km s-i) 


V{5) 
(km s-i) 


J0543+2329 
J0528+2200 


182634 


285 


0.535 


81 


4.435x10"^ 


<5xl0-^" 


0.249 


325±92 


50±131 
51±257 


322±85 
369±178 


J1453-6413 
J1430-6623 


199272 


176 


1.096 


35 


1.756x10-^ 


3.1x10-9 


0.137 


566±147 


-156±257 
319±334 


321±93 
401±118 


J2354+6155 
J2321+6024 


93978 


182 


1.79 


8 


8.513x10-^ 


1.8xl0-y 


0.094 


279±50 


-253±185 
-283±844 


609±102 
930±124 


J1915+1009 
J1909+1102 


222007 


81 


2.646 


10 


4.504x10"'' 


2.02x10-^ 


0.012 


303±64 


-224±306 
-28±59 


616±197 
283±4 


J1832-0827 
J1836-1008 


343000 


55 


3.098 


3 


8.746xl0~6 


1.23x10-^ 


0.00178 


42±17 


20±17 

4±5 


735±47 
2068±630 


J1833-0827 
J1836-1008 


239034 


46 


1.73 


2 


8.367x10-^ 


2.5xl0-y 


0.000753 


26±1 


-360±360 
226±157 


1179±805 
4833±1034 


J1917+1353 
J1926+1648 


221335 


42 


3.399 


2 


9.036x10-6 


<5xl0-9 


0.000593 


696±186 


491±192 
-117±87 


712±346 
402±215 
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p,pc 

Fig.4: Pulsars J0543+2329 and J0528+2200. The 
probability of detecting the pulsar pair in the past 
at a distance of less than p = 5 pc. The solid and 
dashed curves represent, respectively, the probability P 
determined during our simulations and the probability 
Pq of a chance (Poissonian) detection of two (or more) 
pulsars in a cell of space with a size p based on 
the distribution of their density in the Galaxy. The 
nonmonotonic behavior of the function Po{p) at /9>30 
pc is explained by computational noise. The vertical line 
marks po = 285 pc. 

Figure 3 shows an example of such a distribution 
for the pair of J1453-6413 and J1430-6623. Table 
2 presents parameters of the seven pairs: the pulsar 
coordinates, the distance to the pulsar (determined 
from the DM), the proper motion, the characteristic 
age, and the separation between the pulsars at 
the present epoch. For each of them, we then 
simulated 100-300 thousand pairs of trajectories, 
constructed the distributions P{p), Po{p)i and 
T{p), and determined the regions of their closest 
approaches in space (and in the sky), starting from 
p — 5pc. 

For each of the seven pairs, the probabilities 
Po{p) at small p (^50-300 pc) are considerably lower 
than P{p); for very small p, these probabilities differ 
by orders of magnitude, which virtually rules out the 
possibility of a purely chance detection of such close 



encounters in simulations. 

The results of our simulations are presented 
in Table 3. For each pulsar pair, it gives: (1) the 
total number N of simulated trajectory pairs; (2) 
the distance po starting from which the probability 
of a random encounter is higher than that of 
a nonrandom one: Pq(p>pq)>P{p>pq); (3) the 
minimum separation found among the simulated 
trajectories Pmin] (4) the number n of trajectory 
pairs that approach to a distance p ^ 5 pc; (5) the 
corresponding probabilities P (p ^ 5 pc) and Pq 
(p ^ 5 pc); (6) the probability Po(po) of a random 
encounter at the point of intersection between the 
distributions (in the volume of size po); (7) the 
median epoch T and the 50% confidence interval 
that characterizes the approaches to p ^ 5 pc; (8) 
the median values and 50% confidence intervals of 
the current radial velocities Vr for this case; (9) 
the median values of the pulsar peculiar velocities 
V{5) at the epoch of close approach and the 50% 
confidence interval. 



Discussion of Simulation Results 

While analyzing the trajectories of the pulsar pairs 
that may have a common origin, it is interesting (and 
important) to analyze the velocity vectors of their 
separation at CBS disruption, i.e., the magnitudes of 
these velocities and the angles between them, first, 
to test the theories of CBS disruption by a supernova 
explosion and, second, to additionally verify the 
possibility that the pulsars were bound (or unbound) 
in the past. However, the quantities that can be 
obtained in terms of a statistical analysis have their 
peculiarities. Specifically, each individual value of 
the velocity or angle carries little information. In 
our case, these are only individual realizations of 
a random variable that per se do not characterize 
the pulsar pair in any way. Only the characteristics 
of their ensemble are meaningful. In our case, these 
are the confidence intervals, which are wide for these 
quantities. This is understandable: the uncertainty 
in the distances and radial velocities of the sample 
pulsars is very high. Thus, for example, for the closest 
encounters of the pair J0543+2329/J0528+2200 (at 
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/9 ^ 5 pc), the magnitudes of their velocities 
"immediately after the separation" lie with a 50% 
probability in the intervals from ~230 to 400 km 
s~^ and from 190 to 550 km s~^, respectively, while 
the angles between the velocities lie in the interval 
from ~40° to ^^120°. Therefore, it turns out to be 
difficult to reliably estimate the velocity vectors at 
the time of the presumed separation, and such an 
analysis was not performed in this paper. Parallax 
measurements for the pulsars under consideration 
would help greatly to solve this problem. 



J0543+2329 and J0528+2200 

The best results were obtained for the pair of pulsars 
J0543+2329 and J0528+2200 in our simulations of 
182634 trajectory pairs. The approach to 5 pc (or 
smaller) is realized for 81 trajectory pairs, which 
corresponds to the probability P = 4.435x10"'^ 
at a probability of a random encounter Pq for the 
same distances of less than 5x10"^*^. The difference 
by six orders of magnitude suggests that such a 
close encounter of the pulsars is not random and is 
indicative of their common origin. The value of po 
for this pair is 285 pc. For clarity. Fig. 4 shows the 
probability of detecting the pulsar pair in the past 
at a distance smaller than p. (Since the dependence 
is similar in shape for all succeeding pulsar pairs, 
we omit similar figures for them.) The smallest 
separation found between the pulsars in the past is 
0.535 pc. The epoch of close encounters corresponds 
to T ~ 300 000 yr ago, and the characteristic age 
Teh of the younger pulsar J0543+2329 is 253 000 yr, 
which is in satisfactory agreement. If the separation 
between the pulsars is assumed to be at a minimum 
precisely at the birth time of J0543+2329, then its 
actual age is close to 300 000 yr. Since the older 
pulsar J0528+2200 has a characteristic age of about 
1.5 Myr, its progenitor had sufficient time before 
its collapse, in agreement with the evolutionary 
scenario for massive CBSs. The median means of 
the current pulsar radial velocities (for approaches 
to a distance smaller than 5 pc) are 50 and 51 km 
s^^ for J0543+2329 and J0528+2200, respectively. 
The total space velocities of the pulsars relative to 
the LSR at the time of binary disruption (when the 



pulsars approached to p ^ 5 pc) are 322 and 369 km 
s^"*^. Figure 5 shows the pulsar trajectories for which 
the spatial approaches are closest. 



J1453 6413 and J1430 6623 

High probabilities of close approaches are also 
characteristic of the pair of J1453-6413 and J1430- 
6623. Out of the 199272 simulated trajectories, 
35 pairs approached to a distance p ^ 5 pc, 
i.e., the probability P{p ^ 5 pc) is 1.756x10"^ 
and the probability of a random encounter Pq 
is approximately five orders of magnitude lower: 
3.1x10"^. The minimum separation during the 
approach was 1.096 pc (po for this pair is 176 pc) 
(see Fig. 6). The closest approaches occurred about 
570 000 yr ago. The characteristic age of the younger 
pulsar J1453.6413 is about 10^ yr; given the errors, 
this is a good correspondence, as in the previous case. 
The velocities at the binary disruption (when the 
pulsars approached in the past) were about 300 km 
s"^ for J1453-6413 and 400 km s~^ for J1430-6623. 
There is the OB association Pis 20 (Mel'nik and 
Efremov 1995) inside the spatial region of approaches 
to p ^ 10 pc, an argument for the pulsar birth 
and the binary disruption nearby (see, e.g., Bobylev 
2008). The center of the association is located ~80 
pc away from the place in which the approach of 
the pulsars is closest. In this case, the motion of the 
association in the Galaxy is disregarded, since the 
OB association with a velocity much lower than that 
of the pulsar was displaced insignificantly on time 
scales of several hundred thousand years. 



J2354+6155 and J2321+6024 

For the pair of J2354+6155 and J2321+6024 (Fig. 7), 
the minimum separation of the approach is 1.79 pc at 
its probability P = 8.513x10"^. The probability of 
a random encounter Pq is approximately four orders 
of magnitude lower; po for this pair is also fairly 
large: 182 pc. The approaches occurred about 280000 
yr ago, and the characteristic age of the younger 
pulsar J2354+6155 is 920000 yr - this difference is 
also acceptable. 
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Fig. 5: Pulsars J0543+2329 and J0528+2200. Projections of the spatial region in which the measured minimum 
separation between the pulsars does not exceed p = 10 pc. Each point on the plot is the place of closest approach of 
the pulsars for a specific pair of trajectories. The circle and the diamond mark the position of the younger pulsar at 
the present epoch and the position of the older pulsar, respectively. The lines indicate the trajectories at which the 

separation found between the pulsars Pmin is smallest. The asterisk marks the region where it is reached. Here, 

Pmin = 0.535 pc. 
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Fig. 6: Same as Fig. 5 for the pulsars J1453-6413 and J1430-6623. The trajectories at which the separation found 

between the pulsars, pmin = 1.096 pc, is smallest. The points mark the places of the closest approach of specific 

trajectory pairs the minimum separation between which is p ^ 10 pc. 
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Fig. 7: Same as Fig. 5 for the pulsars J2354+6155 and J2321+6024. The trajectories at which the separation found 

between the pulsars, Pmin = 1-79 pc, is indicated by the lines. The points mark the places of the closest approach of 

specific trajectory pairs the minimum separation between which is p ^ 10 pc. 
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Fig. 8: Same as Fig. 5 for the pulsars J1915+1009 and J1909+1102. The lines indicate the trajectories at which the 

separation found between the pulsars, Pmin = 2.646 pc, is smallest. The pointsmark the places of the closest 

approach of specific trajectory pairs the minimum separation between which is p ^ 10 pc. 
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J1915+1009 and J1909+1102 

The results of our trajectory simulations for the 
pair of J1915+1009 and J1909+1102 are presented 
in Fig. 8. Here, the epoch of a close encounter of 
the pulsars corresponds to T '^ 300000 yr. With the 
characteristic age of the younger pulsar J1915+1009 
Teh — 420000 yr, this is a good correspondence. 

J1832 0827 and J1836 1008, J1833 0827 and 
J1836 1008 

In the pairs of J1832-0827 and J1836-1008, 
J1833-0827 and J1836-1008 (Figs. 9 and 10), 
two different secondary components can be 
associated with the same pulsar (J1836-1008). 
Here, we cannot unambiguously choose one of 
these associations based only on our estimates 
of the encounter probability. The two presumed 
companions, J1832-0827 and J1833-0827, are similar 
in characteristic age, both have high velocities (as 
does J1836-1008), and both are at a small distance 
from J1836-1008. Our simulations in both cases 
yielded similar results, and we can only point out 
that both encounters are nonrandom. However, if 
we compare the values of V{5), we will see that 
these velocities are too high for the association 
of J1833-0827 and J1836- 1008. For the pulsar 
J1836-1008, the velocity V{5) reaches almost 5000 
km s~^, a value higher than the velocities of even 
the fastest known pulsars. For the association of 
J1832-0827 and J1836-1008, V{5) for the same 
pulsar decreases to ^^2000 km s^^, a fairly high 
value but not so implausible. Recall that among the 
pulsars with known parallaxes (i.e., with directly 
measured transverse velocities), J1509+5531 has 
the maximum value, which is about 1000 km s^^ 
(Chatterjee et al. 2009). On this basis, it seems that 
the association of J1832-0827 and J1836-1008 is 
more natural. 

However, on the whole, both possible pairs 
suggest a high space velocity of J1836-1008, which 
is marginal even for the pulsar velocity distribution 
in general. This can be an argument against the 
possible association of these pairs. Here, it is hard 



to be guided by the characteristic ages of the pulsars 
(J1836-1008 is the older object in both pairs, and if 
this is the case, then it might not receive a kick at 
binary disruption), because their initial periods are 
unknown. 



J1917+1353 and J1926+1648 

The last pair is J1917+1353 and J1926+1648 (Fig. 
11). The probability P differs from Pq by more 
than three orders of magnitude, suggesting a possible 
association of the pulsars. Their characteristic ages 
are similar, 428000 and 511000 yr, and we cannot 
say which pulsar is actually younger. The encounter 
epoch occurs at T ~ 696000 yr, which may 
correspond to the age of one of the pulsars during 
the formation of which the binary was disrupted. 



CONCLUSIONS 

Out of the initially selected 16 pulsar pairs, 
7 pairs have a common property: the close 
encounters of their components in the past cannot 
be considered random. The selected pulsars have 
enhanced probabilities of close encounters and, 
hence, there is reason to believe that they were 
gravitationally bound in the past, i.e., they were 
members of the same CBS that was subsequently 
disrupted. 

Such disruption is probably the finale of the 
evolution of two initially massive stars — the result 
of the second supernova explosion. When the primary 
(more massive) component during its evolution 
overfills its Roche lobe, mass transfer and outflow 
from the binary begins and it can lose the bulk 
of the hydrogen envelope by the time of collapse. 
As a result, the mass loss by the binary during 
its explosion can be smaller than that needed for 
disruption. In contrast, during the second explosion, 
the binary is disrupted in most cases (Zasov and 
Postnov 2006). 

As was mentioned above, it is virtually 
impossible to identify the case of binary disruption 
after the first explosion, when the secondary 



16 



7 
6 

u 5 

S- 4 

Q 3 










2 




-v • 


1 


1 


1 


1 1 


1 


-8.0 


1 


1 


1 1 


1 


-8.5 


- 




/^ ■ 


- 


-9.0 


- 




f 9 


- 


-9.5 


- 


/ 


„ 


- 


-10.0 


- 


/ 




- 


-10.5 


1 


1 


1 1 


1 



• PSR J 1832-08 27 
^ PSR J 1836- 1008 




279.5 



279.0 



278.5 
a, deg 



278.0 



277.5 



12 3 4 5 6 
D, kfx; 



Fig. 9: Same as Fig. 5 for the pulsars J1832-0827 and J1836-1008. The lines correspond to the trajectories at which 

the separation found between the pulsars, Pmin = 3.098 pc, is smallest. The points mark the places of the closest 

approach of specific trajectory pairs the minimum separation between which is p ^ 10 pc. 
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Fig. 10: Same as Fig. 5 for the pulsars J1833-0827 and J1836-1008. The lines represent the trajectories at which 

the separation found between the pulsars, p„iiri = l-'''3 pc, is smallest. The pointsmark the places of the closest 

approach of specific trajectory pairs the minimum separation between which is p ^ 10 pc. 
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Fig. 11: Same as Fig. 5 for the pulsars J1917+1353 and J1926+1648. The trajectories at which the separation 
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component, being now an isolated one, evolves 
independently from the primary one. Indeed, after 
the CBS disruption through the collapse of the 
primary component, the secondary star will recede 
significantly from the binary position before its 
explosion (in 10^ — 10^ yr). An additional factor 
that breaks the connection between the present- 
day trajectory of the second pulsar and the pair 
coordinates can be an asymmetry in the explosion 
of its progenitor. At the same time, our simulations 
of trajectories is aimed at detecting the encounters 
of pulsars only for a "straight" extension of the 
trajectories into the past. 

Thus, the encounters of pulsars we detected 
probably mark the CBSs that were disrupted during 
the second explosion, while the probability of a 
purely random close encounter between pulsars in 
the past is negligible. 

Note that the differences in characteristic ages 
of the pulsars in the pairs we selected (except for 
the pair of J1917+1353 and J1926+1648) lie within 
the range 5x10^ - 4x10^ yr. These intervals agree 
with the time intervals needed for the evolution 
of the secondary components to be ended with a 
supernova explosion, even if the differences between 
the characteristic and actual ages is taken into 
account. 

We analyzed the relations between the 
characteristic and "kinematic" (determined from the 
epoch of closest approach) pulsar ages. Figure 12 
shows the relationship between these parameters. We 
see no clear correlation; the kinematic age can be 
both older and younger than the characteristic one, 
but the age difference itself is small and comparable 
to the inaccuracy in Teh- 

The technique used here allows us to determine 
whether the kinematic encounters of the selected 
pulsars in the past are a consequence of the natural 
NS density distribution in the Galaxy or such 
configurations are not random. In the former case, 
the pair of pulsars under consideration has never 
constituted a common binary, although the objects 
satisfy the initial selection criteria described above. 
As an illustration of this case, let us consider three 



pulsars: J1835-1106, J1825-0935, and J1824-1945. 
According to the preliminary selection criteria, any 
of the first two pulsars is suitable for being associated 
with the third. However, neither in the pair of J1835- 
1106 and J1824-1945 nor in the pair of J1825-0935 
and J 1824-1 945 can the pulsars be recognized as 
being actually associated. Figures 13 and 14 show 
the probabilities of detecting the components at a 
distance smaller than p. As we see from the figures, po 
for each pair is very small and the difference between 
P and Pq is insignificant. Therefore, in this case, the 
close encounters of the pulsars can be explained by 
purely random factors. 

Thus, we detected six pairs of pulsars that were 
probably members of disrupted CBSs: J0543+2329 
and J0528+2200, J1453-6413 and J1430-6623, 
J2354+6155 and J2321+6024, J1915+1009 and 
J1909+1102, J1832-0827 and J1836-1008, and 
J1917+1353 and J1926+1648. Their further study 
by analyzing a larger number of trajectories will 
allow us to improve our knowledge about the 
pulsars themselves (when the relations between their 
kinematic and "intrinsic" (z/, z>, v) characteristics are 
studied), about the possibility of being members of 
common CBSs in the past, and about the disruption 
kinematics of such binaries. A reliable identification 
of the pulsar pairs will also allow their current radial 
velocities and true ages to be estimated. 
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Fig.13: Pulsars J1835-1106 and J1824-1945. The p dependence of P (solid curve) and Pq (dashed curve). The 

vertical line indicates p^. Here, po is small, while the probabilities P and Pq differ insignificantly. Therefore, the 

conclusion that the pulsars are kinematically associated is invalid. 




Fig.l4: Same as Fig. 13 for the pulsars J1825-0935 and J1824-1945. Just as the pair of J1835-1106 and 

J1824-1945, it is hard to associate the pulsars here. 
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